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ABSTRACT 

A  MACSYMA  program  is  presented  which  determines  whether  a  given  single 
nonlinear  ODE  or  PDE  with  real  polynomial  terms  fulfills  the  necessary 
conditions  for  having  the  Painlevfi  property.  The  program  listing  is  given, 
together  with  a  synopsis  of  the  algorithm,  a  model  of  the  output  of  the 
program,  a  list  of  tested  examples  and  instructions  to  use  the  package. 
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INTRODUCTION 


A  simple  equation  or  system  is  said  to  have  the  Painleve  property,  PP,  if 
its  only  singularities  in  tlie  complex  plane  consist  of  movable  polos'"'^).  This 
requires  that  the  solution  f  of  a  PDE,  say  in  two  independent  variables  (t,x)  can 
be  expressed  as 

OO 

/  =  9“  E  (*) 


•iSPECTEO  . 

^  *  y 


with  uo{t,x)  7^  0,  a  a  negative  integer,  and  where  tifc(Ga:)  are  analytic  functions 
in  a  neighbourhood  of  the  singular,  non-characteristic  manifold  g{t,x)  =  0,  with 
9x{t,x)  #  0. 

Performing  the  Painleve  test,  as  proposed  by  VVeiss^\  involves  three  steps: 

(i)  Determination  of  the  negative  integer  q  and  Un  from  the  leading  order  ansatz 
f  a  uo?“,  by  balancing  the  minimal  power  terms; 

(ii)  Indentification  of  the  non-negative  integer  powers  r,  the  resonances,  at  which 
arbitrary  functions  Ur  can  enter  the  expansion  /  oc  uo<7“  4  Urff"'*"’’.  This  is 
achieved  by  requiring  that  Ur  is  arbitrary  after  substitution  of  this  form  for  f  into 
the  equation,  only  retaining  the  most  singular  terms; 

(iii)  Verification  that  tin*  correct  number  of  arbitrary  functions  Ur  indeed  exists 
by  substituting  a  truncated  expansion  of  the  form  (I),  k  =  1,2, ,..,rmax,  where 
rmax  represents  the  largest  resonance,  into  the  given  equaton.  At  non-resonance 
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levels,  Ufc  should  he  uniUiili'iKuousIy  delerniined;  at  resonance  levels  should  be 
arbitrary  due  to  a  vanisliing  coeHiciciit  of  (compatibility  condition). 

Here,  minpowg  denotes  tin*  (negative)  power  in  g  of  the  most  singular  terms  in 
the  equation. 

An  equation  or  system  for  which  the  above  steps  can  be  carried  out  consis¬ 
tently,  is  said  to  have  the  IM*  and  is  conjectured  to  be  inlegrable'''’ '  ('c  inter 

examples^ “ of  integrable  I’DHs  without  the  IM’  disprove  the  necessity  of  the 
PP  for  integrability  (i.  e.  lu'ing  exactly  solvable  by  a  linear  integral  equation  of 
Inverse  Scattering  Transform);  whereas  some  non  -  integrable  equations  seem  to 
have  the  PP  as  defined  by  Wei.s.s  et  al^\  hence  questioning  the  conjecture'*^  in 
its  present  form.  Refinements  of  the  PP  have  been  established*^”*®^  allowing  for 
rational  power  expansions  of  f,  which  lead  to  algebraic  branch  points  in  addition 
to  movable  poles. 

Nevertheless,  the  PP  has  been  a  useful  criterion  for  complete  integrability 
tests®’®’*°”'®^.  It  is  a  handy  tool  in  the  derivation  of  Hacklund  transformations 
and  Lax  pair  representations  for  various  PDEs®  *®“^°\  and  it  gives  insight  in 
Hirota’s  direct  method  for  solitary  wave  construction^* and  other  expansion 
methods  which  lead  to  rational  solutions*®^. 

2.  SCOPE  AND  LIMITATIONS  OF  THE  PROGRAM 

The  program  is  largely  based  tipon  the  M.ACSYMA  routine  ODKPA INLEVE 
developed  by  Rand  and  Winternitz^'^^  which  checks  the  PP  for  a  single  ODE  with 
real  polynomial  terms. 

The  package  works  for  both  single  ODEs  and  PDEs,  with  arbitrary  degree 
of  nonlinearity.  In  principal,  the  number  of  variables  is  unlimited.  Although  the 
released  version  of  the  program  works  with  only  four  independent  variables.  PDEs 
may  have  time  and  spare  ilependent  coefficients  of  integer  degree.  'I'ransendental 
terms  in  the  equation  must  be  removed  by  a  suitable  exponential  transformation 
of  the  dependent  variable’^‘^‘’'^®L 

Due  to  ’intermediate  e.xpre.ssion  swell’  it  jiiay  occur  that  the  program  can 
not  perform  the  test  in  a  reasonable  time.  Therefore,  verifying  the  resonances 
and  compatibility  condilioiis  are  made  optional  through  the  boolian  variable 
do -resonances.  The  higlu'st  level  of  verification  is  controllable  by  entering  a  value 
for  max.resonance.  For  the  investigation  of  Macklund  transformations  one  may 
consider  calculations  beyond  the  automatically  determined  level  rtnax  by  taking 
max-resonance  >  rrnai. 

Intermediate  output  ran  be  obtained  by  including  extra  print  statements  in 
the  program.  (Warning  ;  After  any  alternation  in  the  main  program  a  new  LISP 
version  needs  to  be  created  by  bat ch (main. p)  in  .MACSYMA). 


Major  dinicultirs  aiisc'  wlu'ii  u^,  can  not  be  substituted  since  it  occurs  in 
subsequent  calculations  at  a  lower  def^ree  than  the  one  present  at  its  evaluation. 
In  that  case,  calculations  l»y  liand  or  in  intiTactive  M  ACISYMA  are  reconmiended. 
The  vital  steps  for  tlial  aie  easily  e.xtracted  from  the  program  listing. 

At  present,  the  program  does  not  perform  the  weakest  forms  of  the  I’ainleve 
test*^“'^^  although  a  in  (I)  may  l)0  rational  and  can  be  positive,  through  a  user 
supplied  value  for  beta.  I^atcr  extonsions^^’^  will  cover  complex  equations  and 
coupled  systems.  Some  of  these  must  be  treated  by  rather  general  rational  power 
expansions’*”*®^.  A  further  re/inemenl  of  ihe  algorithm,  performing  selective 
substitutions  of  Uk  and  their  derivatives  in  the  recursive  relations,  will  allow  to 
construct  Backlund  transformations  for  PDEs®’’®’^®\ 

3.  HOW  TO  USE  THE  PROGRAM  :  EXAMJ’LES 

The  equation  to  be  tested  should  be  saved  in  a  batch  file,  before  entering 
MACSYMA. 

Example  1  :  An  ODE,  where  the  use  of  independent  variable  x  is  mandatory. 

/*  Batch  file  R.W.(5.G)  lor  testing  Eq.  (5.6)  in  paper  Rand-Winternitz^^^  */ 
eq;  (f '2)*fxl3l(x)  -  3^(fxill(x))*3  $ 
beta  :  -2  $ 
max_resonance  ;  5  $ 

/*  beta  and  max-rcsonaiice  arc  specified  here  */ 

Note  that,  during  the  lest  the  variable  x  will  be  replaced  hy  g  —  x  -  iq)  where 
xq  is  the  arbitrary  initial  value  of  x. 

Example  2  :  A  PUE,  with  independent  variables  t,x,  and  y. 

/*  Batch  file  KP,  for  the  Kadomtsev-Petviashvili  Equation^^  */ 

eq  :  ftxy|l,l,0|(t,x,y)  +  (ftxy|0,l ,Ol(t,x.y))  ‘ 2  f  P ftxylt),2,()j(t,x,y) 

+  b*ftxy(0,4,0)(t,x,y)  I  ftxy|0,0,2){l,x,y)  S 
/*  b  is  a  constant;  beta,  do  resonances  and  max  resonance  are  not  specified  * / 

Similarly,  using  // j-ycj., .|(t,  x,  y,  z),  evolution  and  wave  ecpiations  in  t,  x,y 
and  z  can  be  entered,  for  pararru’trical  equations  the  explicit  dependence  on  the 
variables  must  be  given  (e.g.  coefficients  o{t),  6(x), ...).  The  five  first  letters 
of  the  alphabet  are  reserved  to  denote  arbitrary  constants,  if  necessary  to  be 
supplemented  with  a  1,  u2, ...,  Al , .... 


Once  the  batch  files  are  marie,  enter  MAC'SVMA  and  subsorpiently  type  : 


(cl)  batch(main-p)  $ 
(c2)  batch(selup-p)  $ 
{c3)  batch(kp)  $ 

(c5)  writefile(test-kp)  $ 
(c6)  batch(exec_p)  $ 
(c7)  closefile()  S 

4.  TEST  RUNS 


to  create  pamleve.lsp,  only  to  be  doiu'  ones 
*  to  repeat  before  every  new  exaiiiple 


— ►  reads  in  the  batch  file  kp 
-*  opens  output  file  test  kp 
— »  starts  the  actual  I’ainleve  test 
— »  closes  output  file  test-kp 


Although  the  program  has  been  successfully  tested,  with  eunice  MACSYMA 
version  300.3  on  a  VAX  1 1 /7H(),  for  more  than  30  (>xamples  including  all  ODEs  and 
selected  PDEs  from  the  citerl  references,  uircxpected  situations  may  still  occur. 
In  case  of  trouble,  the  printout  gives  vital  information  to  remedy  the  problems. 
The  program  does  not  rely  on  non-standard  MAC'SYMA  routines  (e.g.  powers  is 
provided). 


Model  of  Output  : 

For  the  Space  Dependent  llurgers  equation :  -x^/2x  -  i/Zi  d  ft  '=  0,  which 
passes  the  test  : 


2 

PAINLEVE  ANALYSIS  OF  EQUATION.  -  f  *  -  f  f  *  +  f 

K  X  X 


alia 

SUBSTITUTE  u  g  FOB  f  IN  ORIGINAL  EQUATION. 

0 

MINIMUM  POWERS  OF  g  ARE  [2  alfo  -  1,  olfo  -  2] 


e 


2  alfo  -  1  2 

•  COEFFICIENT  OF  g  IS  -  u  olfo  g  x 

0  X 

olfo  -  2  22 

«  COEFFICIENT  OF  g  IS  -  u  (olfo  -  1)  olfo  (g  )  x 

0  X 

FOR  EXPONENTS  (  2  olfo  -  f  )  AND  (  olfo  -  2  )  OF  g,  DO 
WITH  olfo  -  -  1  .  POWER  OF  g  is  -  3  — >  SOLVE  FOR  u 


TERM  -u  g  x(2g  x-u)—  IS  DOMINANT  IN  EQUATION. 
0  X  X  0  3 


9 

WITH  u  -  2  g  X  — >  FIND  RESONANCES 
0  X 

olfo  r  +  ol f o 

SUBSTITUTE  u  g  +  u  g  FOR  f  IN  EQUATION 

0  r 


^  2  2  r  -  3  . 

TERM  -  (g  )  (r  -  2)  (r  +  t)  x  u  g  IS  DOMINANT  IN  EQUATION. 

*  r 

THE  ONLY  NON-NEGATIVE  INTEGRAL  ROOT  IS  [r  -  2] 

WITH  MAXIMUM  RESONANCE  -  2  — >  CHECK  RESONANCES. 


\  k  -  1 

SUBSTITUTE  POWER  SERIES  >  g  u  FOR  f  IN  EQUATION 

/  k 

k  *  0 

WITH  u  -  2  g  X 
e  X 

1  2 

COEFFICIENT  OF—  IS  2g  x(g  x  +u  g  x-g) 

2  XXX  lx  t 

9 


g  K  -  g 

X  X  t 


COEFFICIENT  OF  -  IS  0 
9 

u  IS  ARBITRARY  I 

2 

COMPATIBILITY  CONDITION  IS  SATISFIED  1 


The  results  of  test  runs  on  15  other  examples  are  summarised  in  Table  1. 

5.  LISTING 

5.1  Listing  of  setup.p  routine 

/•  The  botch  file  SETUP_P,  which  initializes  the  program  after  every  example  •/ 
ki I  I  (al  l)$ 
loodprint  :  falseS 
loadfile  ("painleve.  Isp'')$ 
derivabbrev  :  true$ 
depended  f. eg], It, x,y,z))$ 
fx[k](xj:»diff(f.x,k)$ 
ftx[k, l](t,x):-dlff(diff(f ,t.k).x.  I)$ 
ftxy[k. I .mi(t .x.y) ;-di f f (dif ((diff (f . t .k) ,x. I) .y.m)$ 
ftxyz[k,l,m,ni(l,x,y,z):-dlff(diff(diff(d:ff(f.t,k),x,l),y,m).z.n)$ 


5.2  Listing  of  tirr  p  routine 


/•  The  batch  file  EXEC_P,  to  execute  the  test  for  an  example  •/ 
exec_palnleve  (eq,  beta.  do_resonances ,  max_resonance)$ 


5.3  Listing  of  matn.p  routine 


/•  Tha  botch  file  MAIN_P,  with  the  moin  program  for  the  Poinleve  test  •/ 
def i ne_vor I  obi e  (do_resonance3.  true,  boolean)$ 
declare  (max_resonance ,  lnteger)$ 

nonposint  (expres)  ;«  If  Inf egerp(expros)  then  I s(expres<-0)  else  false$ 

so  I veforh I pow  (expr.vor)  block(  [m], 
expr  ;  expand(expr  1  , 
m  ;  h I pow (expr.vor), 
return(port(3olve(expr.vartm),1)))$ 

konstont  (expr.vor.q)  block  ((h,co), 

CO  :  rot coe f (expond( expr ), vor , 0) , 

I f  q  I  1  then  ( 

CO  ;  rotsubst  (h . vor t ( 1/q) ,co) , 

CO  ;  rotcoef {co,h.0) 

). 

return  (cc))$ 
pow(exp r , vor ) ; «bl nek  ( 

(power_l I  8 1 :  ‘ [ | . f  f : expand (expr) , tmp, tmpi , Infun, Inf  log: true], 
dec lare(power_l ist .spoclol), 

if  not  otom(f7)  ond  part(ff,0)  "+"  then  I nf un ; orgst f f )  else  lnfun:[ff], 
for  i i^  In  infun  do  ( 
tmpi : I iX, 

if  not  atom(lmpl)  and  port(tmpl.0)  •  then  tmpi : substpor t ("+”, tmpi ,0) , 

it  otom(tmpl)  or  inpart(tmp1 ,0)#'’»“  then  tmp1:[tmp1], 
tmp .mop  I i st ( I umbdo( lu] ,qet_exponent (u, vor ) ) , tmpi ) , 
tmp : de I ete( ' [ ] , tmp ) , 
i f  tmp  •  ' [ J  t  hen  ( 

If  not  member(0,power_l let)  then  power_l ist :endcons(0,power_l ist) 

) 

else  ( 

if  not  member ( f I  re t ( tmp) , power_ I  I s t )  then 
power.,  I ist: endcone( f I rst ( tmp) , power. I ist) 

) 

). 

power. I i st )$ 

get.exponen t (expr , vor ) :"bl ock  ( 

if  f reeof ( VO r , expr )  then  retorn(’[]), 
if  otom(expr)  then  return(l), 

if  inport(expr,0)»"t"  ond  i nport (expr , 1 )*vor  then 
return( inpart(expr,2)) 
el  se 

return(’(]))$ 

check.rosononce3(eq , VO  1 0  I f o ,ml npowg , uzero, rmox , s)  :■  block  ( 

[series, eql ,eq2,dnl fa.co,uk), 

print  ("  WIIH  MAXIMUM  RESONANCE  -  ",  rmox,  "  — >  CHECK  RESONANCES."), 
print  ("  SUBSri lUIE  POWER  SERIES  ",  ’ 8um( u[ k] egT ( vo I o I f o+k) , k . 0 , rmox) . 

"  FOR  f  IN  EQUATION."), 
series  :  sum( u[ k | xg t ( vo I o I to+k) , k , 0 , rmox) , 
eql  ro t subs  I ( ser I ee , f , eq) , 
eql  :  ev( eql , d I f  t ) , 
eql  :  expand(eql  /  ( g tml npowg) ) . 
eql  ro t  s  i mp(  ro  t  subs  t  (0,  gt  ( s* rmox-kl  ), eql  )) , 

CO  :  denom(eql), 

if  Integerpfco)  ond  co  |  1  then  eql  :  num(eql), 

if  uzero  f  []  then  eql  :  rot subst ( rhs(expand(uzero) ) . I hs(uzero) , eql ) , 


6 


iolfo  :  denom(oba(volal (a)K 
:o  ;  konstant  (eq1 ,g,dal fa) , 
if  CO  I  0  then  ( 

print ("•••  FATAL  ERROR  •  NON-ZERO  LEADING  COEFFICIENT  REMAINS  ", 

"AFTER  SUBSTITUTING  ",  uzoro.  ".  •••"), 
return  () 

it  dalfo  I  1  then  eql  :  radcan(eql), 
if  uzero  f  []  then  print!"  WITH",  uzero) 

else  prlntr*  WITH".  u[0].  "ARBITRARY  I"), 
for  i  1  thru  s*rma*  do  block  ([k], 

if  dolfo  I  1  then  (or  j  :  T  thru  dolfo-l  do  ( 

CO  ;  rotcoef (expnndCeql ) , q, j/dol to) , 
if  CO  I  0  then  ( 

if  not  freeo((q,co)  then  ( 

printC'***  FATAL  ERROR  :  COEFFICIENT  STILL  CONTAINS  g.  •••"). 
return 

) ,  /•  end  i (  •/ 

eql  :  radcan(eq1  -  co*gt ( j/do I f a) ) 

)  ). 

eql  ;  ratsinip(eq1/g)  . 

CO  :  konstont  ( eql . g , do  I f a) , 
eql  ;  eq1  -  co. 

CO  :  rots imp( xthru(co) ) . 

if  not  integerp(co)  then  co  :  factor(co). 

print("  COEFFICIENT  OF”  ,  gf  (  l-Hit  i  npowg) ,  ”  IS  ",co), 

k  :  i/s, 

i f  CO  #  0  and  integ«rp(k)  then  ( 
if  rotcoef  (co ,  u[  k  )  .  s)  ••  0  then  ( 
if  not  ( reeof ( u[ k J . co)  then 

print("»««  RFSONANCE  VIOLATED  BY  APPEARANCE  OF".u[k], 

"FROM  NON-LEAOINO  TERMS.  •••") 
else  b(ock(  [top, J ,cox0l, 

print("  ",  u[k)t8.  "  IS  ARBITRARY  I"). 

print^•  COMPATIBILITY  CONDITION:  ",  co.  "  -  0"). 

if  (rotsimp(ev(co.  diff))  -  0)  then 

print("  CONDITION  IS  SATISFIED  I") 

else  ( 

print!"  •••  CONDITION  NOT  SATISFIED.  •••"), 
print("  •••  CHECK  FOR  FREE  PARAMETERS  OR  PRESENCE  OF". 
u[0l,".  •••")), 
top  :  num(co) , 

if  not  f reeof (*0. top)  then  for  j :0  while  top  |  0  do  ( 
top  :  expond(top), 
co*0  :  coef f ( t op . *0.0) , 
top  :  ( top-co*0)/*0, 

COK0  :  factor(coxO) , 

printC  COEFFICIENT  OF  ".x0tj,"  IN  NUMERATOR  IS".cox0) 

)  )  )  /•  snd  do  «/ 

else  block  (rtempl. 

uk ;  solvefco.ii  kits), 
uk  :  part(uk, I ) , 

temp  ;  f oc t or ( ev ( rha( uk) .dl f f )) , 
uk  :  rot  subs t ( t emp , rha(uk) , uk) , 
pr i n t ( "  "  ,  uk ) . 

eql  ra t subs t ( rhs( uk ) , I he(uk) , eql ) , 

if  uzero  |  []  then  eqt  ;  rot  subs t ( rhs(expand(uzero) ) , I h8(uzero) ,eq1 ) , 
eql  :  ev(eq1 ,dl f  f ) 

)  )  /•  end  else  If  •/ 

else  if  CO  ••  0  and  Integerp(k)  then  ( 


print  f"  ",  ulkJtB.  "IS  ARBITRARY  I"), 

print  ("  COMPATIBILITY  CONDITION  IS  SATISFIED  I") 

)  ))$  /•  end  of  function  chock_resononce9  •/ 


f  i  nd_re90nance9(eq  .  VO  I  a  I  f  a  ,m  1  npo«<j  ,  UT«  ro)  block  ( 

[eql ,eq2,coef  r,9. I . I  I . tempt ,t  emp2 . powUB , do  I f o  ] , 

print  ("  WITH  ".uzero,  "  — >  FIND  RESONANCES"). 

print  ("  SUBSTITUTF  ",  u[0]«gfalto.  "  +  ",  u[ r ] , gt (a I f o+r) . 

"  FOR  f  IN  EQUATION"), 

eql  ro 1 9ub9 1 (u[0) • ( g t VO  I n I f a)+u[ r ] • (gt ( r+vQ I o 1 f o) ) , f , eq) . 

eql  :  ev(eq1 , d i f  f )  , 

eql  :  rat9imp(eq1/(gfml npoog)) , 

dalfo  :  denom  (ab9( vo I o I f o) ) , 

if  dalfo  I  1  then  eql  rndcon(eql). 

powU0  :  h  i  pow(  I  h9  ( uio  ro)  ,  11(0  ] )  , 

tempi  :  denom(eql),  if  I n I ege r p( t empi )  ond  tempi  |  1  then  eql  :  num(eql), 
temp2  :  h i pow( eube t ( r-1 , eq 1 ) . g) . 
coefr  :  0. 

for  i  ;  1  thru  temp2  wiiile  coefr  -  0  do  ( 
eq2  ;  rotcoef (eql . u( r  ]  ,  I  )  . 
eq2  :  rotcoef (eq2 , g . i • r ) , 
coefr  :  konetont  (e<|2 ,  g ,  do  I  f  o)  , 
i f  uzero  #  [ ]  then  ( 

coefr  :  9xpand(coe f r )  , 

while  h i pow( coe f r , u( 0 ) )  >-  powU0  do 


coefr  :  expand( rot eubet ( rh9(eKpand(uzero) ) ,  lh9(uzero),  coefr)) 


coefr  ;  f oc t or ( num( ro 1 8  I mp( rodcon(coef r/u[ r ] 1 9) ) ) ) , 

print  ("  TERM  " , coef r . u( r J t 9 , gt ( 9* r+m i npowg ) , "  IS  DOMINANT  IN  EQUATION."), 

coef  r:xthru(coefr/factor(qcd(coofr,diff(coefr,r),r))), 
if  port(coef r,0)  »  then  coefr  :  -  coefr, 

if  port (coef r .0)  *  "/"  then  coefr  :  num(coefr), 
temp2  :  1 . 

if  port (coef r .0)  -  "+"  then  ( 
if  hipow(coef r. r)  =  1  then  ( 

if  (coef f (coef r . r)  ■  1  and  nonpo9int(9ub9t(r>0,coef r)}}  then 
temp2  :  coefr 

)  ) 

elee  if  port (coef r .0)  -  then 

for  ii  :  1  thru  I  eng t h( coe f r )  do  block(  [pp], 
pp  :  port (coef  r , I  i ) , 
if  hipow(pp.r)  -  1  then  ( 

if  coeff(pp,r)  -  1  and  nonpo9 I nt (eube t ( r»0, pp) )  then 
temp2  ;  temp2  *  pp 

)  ). 

I  f  temp2  «  1  then 

print  (”  THERE  ARE  NO  NOM-NEGATIVE  INTEGRAL  ROOTS  FOR  r.") 

elee  ( 

temp2  :  eol ve( temp2, r) . 
if  Iength(temp2)  =  1  then 

print  ("  THE  ONLY  NON-NEGATIVE  INTEGRAL  ROOT  IS  ",  temp2) 

e  I  ee 

print  ("  THE". Ienqth(temp2). "NON-NEGATIVE  INTEGRAL  ROOTS  ARE  ",  temp2). 

print  ("  "). 

temp2  :  opp I y (max ,map( rhe , temp2)) . 
if  temp2  >  0  then  ( 

if  do_re9onance9  then  ( 

it  i ntegerp(mox_re9onance)  then  temp2  :  max_re9ononce , 


chack_re9onance*(aq , va I  a  I f a .mi npoarg . uz«ro , temp? , s) 

) 

el  ee 

print  ("  EXIT  E INO.RESONANCES  WITHOUT  CHECKING.") 

)  ))$  /•  end  of  function  f I nd_re9onances  •/ 

f i nd_uzero(eq . VO  I  0  I f 0 .m I npowq)  -  block(  [eq? , r , I  owes t , uzero . do  I f a] . 
print  (■'  WITH  olfo  -  ",  votolfo,  ",  POWER  OF  g  is  ",  minpowg, 

"  — >  SOLVE  FOR  ul0]K 
eq?  :  rot  subs t ( u[0) • ( q f vo I o I f a ) , f , eq) , 
eq?  :  ev(eq2 . d  i  f  f ) . 

eq?  :  rotsimp(eq2  /  ( q tm I npoeg) )  .  ' 

dolfo  :  denom  (abs( vo I o I f o ) ) , 
if  dolfo  I  1  then  eq?  :  radcon(eq?), 
lowest  :  denom(eq?). 

if  i nteqerp( I owest )  and  lowest  |  1  then  eq?  :  num(eq?), 
lowest  f oc t or ( konst  on t  (eq? , q , do  I f o) ) . 

print  ("  TERM  ".lowest,  qfminpowq.  "  IS  DOMINANT  IN  EQUATION."), 

if  lowest  =  0  then  f I nd_reeononces  ( eq . vo I o I f o .m i npowg , [ ] ) 
else  if  f reeof (u[0] , I owest )  then 

print("*«»  FAItllRF  OF  PAINLEVE  TEST  —  NONZERO  LEADINGORDER  IS  ", 
"FREE  OF".u[0].".  ABANDON  THIS  olfo  VALUE.  •••") 

else  ( 

for  i:1  while  remo i nde r ( I  owes t ,u(0] )  ”  0  do  lowest  :  x t h r u( I  owes t/u[0] ) , 
if  f reeof ( u[0] . I  owes t )  then 

print("  •••  FAIIURE  OF  PAINLEVE  TEST.  ".  u[0]-0. 

".  ABANDON  THIS  olfo  VALUE.  •••") 

else  ( 

if  (denom( I  owes t )  §  1)  then  lowest  :  num(lowe9t). 

lowest  :  f oc tor ( I  owes t/(qcd( I owest . d i f t ( I  owes t . u[0) ) , g[0) ) ) ) , 

if  port( lowest .0)  -  "-"  then  lowest  ;  -lowest. 

if  par t ( I  owes t . 0)  «  "/"  then  lowest  ;  num( lowest), 

if  part( lowest ,0)  f  then  ( 

uzero  :  sol vef orh i pow( I owest . u(0] ) , 

uzero  :  rotsub9t(faetor(rh9('jzero)).rhs(uzero), uzero) , 
f ind_re9onance8(eq,valal fo.mi npowg, uzero) 

) 

else  for  1:1  thru  length( lowest)  do  ( 
uzero  ;  par t ( I  owes! . i ) . 

If  not  f reeof (u[ . uze ro)  then  ( 

uzero  ;  Bolveforhipow(uzero.u[0]) , 
flod_re9onance9(eq.valolfa.tiii  npowg  .uzero) 

)  )  )  ))$  /«  end  of  function  find_uzero  •/ 

poinleve  (eq.  beta.  do_r eeononces ,  mox_re9ononce)  block  ( 

[eqt,  u0,  olfo,  I,  exp,  list,  len,  eq?,  olfolist], 

print  (" - - - — 

print  ("PAINLEVE  ANALYSIS  OF  EQUATION,  ",  eq.  "  -  0"). 

print  (" - 

print  ("  SUBSTITUTE  ".  u(0]*qtalfa.  "  FOR  f  IN  ORIGINAL  EQUATION."  ). 

eqt  :  rot  subs t ( ul 0] vq t a  I f o , f , eq)  , 

eqt  :  ev( eq 1 . d i f  f ) . 

eqt  :  fiiap(  f  oc t or  .  eq  1 )  . 

i  ;  h i pow(expand(denom(eq1 ) ) , g) , 

list  :  pow(num(eq1 ) ,q) . 

list  ;  map( I ambda( [ent ry ] .en t ry-l ) , I  I st ) , 
list  :  opp I y (m i n . I i st  )  , 

if  part(tist,0)  <•  min  then  list  ;  args(list)  else  list  ;  [list], 
len  :  length( list). 

print("  MINIMUM  POWERS  OF  g  ARE  ".list). 


0xp  :  expond(eql), 
for  i : 1  thru  Ian  do  ( 

eq2  :  expand( ro t coef («xp , g , par l( I  I  at , I ) ) ) , 

aq2  :  rotcoef {eq2 , g ,0) , 

eq2  :  rat  subs t (0 , gta I ( a , aq2) , 

eq2  ;  factor(eq2) . 


print(‘ 


•  COEFFICIENT  OF  » . gtpar t ( I i st , i ) .  "  IS  ",  eq2) 


ol f al i st : [  ] , 

if  Ian  >  1  and  not  rotnump  (bata)  then 
for  i  1  thru  len-1  do  for  j;i+1  thru  Ian  do 
blocl<(  [  vo  I  o  I  f  a  .  I  i  s  1 1  .m  I  npowg  ]  . 

print  ("  — - "). 

vololfo  :  so  I  ve( po r t ( I  I s t , I )-oa r t ( I  I  a t , J )  ,  a  I  f  o)  , 
volalfo  :  rhs( pa r t ( VO  I  a  I f a , I ) ) , 
if  volalfo  <  0  t  han  ( 

listi  ;  subs t ( o I f o-vo I o I f o , I  I  a t ) , 

Hiinpowg  :  part(llaM,l), 

if  not  friamber(volal  fn.nl  fal  let)  than  ( 

print  ("  FOR  EXPONENTS  ( "  .  por  t  (  I  i  s  t .  I  )  .  *' )  AND  ("  .  port  ( I  i  at .  j )  ,  " )  OF  g,  DO") 
olfalist  :  cons ( VO ! a  I f 0 . o I f o I i a t )  . 
if  app I y (m i n . I i st I )  |  mlnpowg  then 

print("  POWER  OF  g  IS  NOT  MINIMAL  —  SKIP  THIS  VALUE  OF  ALFA") 
else  f ind_ozaro(aq.volal fo.mlnpowg) 

) 

el  se 

print  ("  FOR  EXPONENTS  ( "  .  par  t  (  I  i  st .  I ) .  *' )  AND  (  "  .  por  t  (  I  i  s  t ,  j  )  . "  )  OF  g.". 

"  alfo  -  volalfo,  "  ALREADY  DONEI") 

else  print("  FOR  EXPONENTS  (" . par t ( M st . i ) . ")  AND  ( ” . por t ( I  I s t , j ) . 

")  OF  g.  olfa  -  ".volalfa."  IS  NON-NEGATIVE.") 

),  /•  This  is  the  end  of  the  first  line  including  the  block  •/ 


list  ;  subs t (be to . o I f a , I  I st ) , 
list  :  opp I y (m i n , I i s t )  . 
if  ratn'jmp(bata)  then  ( 

print("  RESULTS  FOR  USER-SUPPLIED  VALUE  olfo  -".beta), 
f i nd_uzero(eq , beto , I  let) 

) 

)$  /•  end  of  function  palnlove  •/ 

exac_pa i n I  eve  (eq.  beta,  do_reaononca8.  max_rosononco)  block  ( 

[degpdf  t , degpdf  x . degpdf  y , degpdf  z] , 
degpdft  :  dar i vdegrae(aq , f  ,  t )  , 
degpdfx  :  do r i vdeg rae( eq , f , x )  , 
degpdfy  :  da r i vdeg ree(aq . f . y ) , 
degpdfz  :  dar i vdegree(eq , f . z ) , 

if  (degpdf t»=0)  and  (dagpdfy-0)  ond  (dogpdfz»0)  than  block  ( [ I , degx ,aql ] , 
romove( [ f , g . eq] , dependency ) . 
depends( f .g) . 
dopends(g. x) , 
depends(u , k)  , 
declaro(g, mo  invar), 
aql  :  ev ( eq . d I f  f ) , 
eql  ;  rotsubst(l,diff(g,x),oqt), 
degx  :  der I vdegree(eq I , g , x) , 
for  i ; 2  thru  degx  do  ( 

aql  ;  rot  subs t (0 , d I f f (g , x ,av( I ) ) ,aq1 ) . 
eql  :  num(xthru(eql ) ) 


eq  :  subst(g4x0, x ,eq1 ) , 

print  ("  SUBSTITUTE  x.  g+x0) 

) 

•  I  se  ( 

<Jepen<^->  ([a].tt.x,y.tl). 
depenq9( [uj  .[k.t.x.y.i]) 

painl«v«  (eq,  beto,  do_reeonanc«9 ,  max.resononce) )$ 
*ove  ( ( "po i nieve . I »p"] ,  do_re«ononce»,  functions); 
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